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VFL/LF  Reflection  Properties  of 
the  Low  Latitude  Ionosphere 


1.  INTRODUCTION 

This  report  summarizes  data  describing  the  longwave  reflection  properties  of  the  iOw-latitude 
ionosphere  in  Brazil,  obtained  during  an  experimental  program  conducted  jointly  by  the  Air  Forces  of 
the  United  States  and  Brazil.  The  data  were  obtained  during  May  and  June  of  1980  in  southern  Brazil 
using  a  unique  VLF/LF  pulse-ionosounding  system  operated  by  the  Brazilian  Air  Force.  The 
transmitter  of  the  ionosounder  was  located  near  Paula  Freitas  (see  Figures  1  and  2).  and  three  portable 
receiving  systems  were  used  to  obtain  measurements  at  nine  locations  to  the  east  and  west  of  it.  at 
distances  ranging  from  156  km  to  280  km,  as  indicated  in  Figure  2. 

2.  VLF/LF  IONOSOUNDING  TECHNIQUE 

Details  of  the  VLF/LF  pulse-ionosounding  technique  are  given  in  a  number  of  reports  describing 
its  use  in  other  experiments  conducted  at  low-,  mid-,  and  high-latitudes  (Rasmussen  et  al1.  Lewis  et 
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Figure  1.  Brazilian  Air  Force  High  Resolution  VLF/LF  lonosounding 
Transmitter  Site  At  Paula  Freitas,  Brazil. 
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al2,  and  Kossey  et  al3,  respectively).  Figure  3  Illustrates  the  technique.  The  transmissions  consist  of  a 
series  of  extremely  short  pulses  (55  nsec  for  the  experiment  in  Brazil),  precisely  controlled  by  a  cesium 
beam  time  standard,  and  radiated  from  a  vertical  antenna  (130  m  in  Brazil).  The  transmitted  pulses 
are  so  short  that,  for  distances  out  to  several  hundred  kilometers,  the  groundwaves  arrive  and  are 
recorded  before  the  arrival  of  the  ionospherically  reflected  skywaves.  Thus,  data  are  received  free  of 
ambiguities  arising  from  groundwave  and  skywave  interference.  The  transmitted  pulses  are 
vertically  polarized,  but  because  of  the  effects  of  the  geomagnetic  field  the  reflected  skywaves  are 
iliptically  polarized.  Separate  loop  antennas  are  used  to  sense  the  so-called  "normal"  (II  -)  and 
"converted"  (.L-)  components  of  the  downcoming  skywaves.  One  antenna,  with  its  axis  oriented  in  the 
plane  of  propagation,  senses  the  groundwaves  and  the  normal  components  of  the  skywaves,  while  the 
other  antenna,  orthogonal  to  it,  is  nulled  on  the  groundwaves  and  senses  the  converted  components  of 
the  skywaves.  Examples  of  the  loop  antennas  used  in  Brazil  are  shown  in  Figure  4. 

Figure  5  gives  typical  examples  of  groundwave  and  skywave  pulses  received  in  Brazil.  The 
separation  in  time  between  the  normal  skywave  pulse  and  the  groundwave  pulse,  which  is  a  feature  of 
the  ionosounding  technique,  is  illustrated  clearly  in  the  figure.  Figure  6  shows  the  Fourier  spectrum 
of  a  transmitted  pulse  to  illustrate  the  wide  range  of  VLF /LF  frequencies  which  were  available  to  probe 
the  ionosphere.  In  practice,  the  pulses  were  transmitted  at  a  precisely  controlled  repetition  rate  of  399 
pulses/sec,  so  that  the  spectrum  of  the  received  pulses  consisted  of  discrete  frequency  components,  399 
Hz  apart. 

As  described  by  Lewis  et  al2  (op.  cit.),  after  Fourier  analyses  of  the  received  pulses,  the  group 
delays  between  the  skywaves  and  groundwaves  can  be  used  to  determine  effective  heights  of  the 
ionosphere;  and,  the  relative  amplitudes  of  the  skywave  and  groundwave  frequency  components  can 
be  used  to  determine  magnitudes  of  the  normal  (||R  h)  and  converted  (i|Rj_)  reflection  coefficients  of  the 
ionosphere4.  The  experiments  In  Brazil  yielded  ionospheric  reflection  heights  and  effective  plane- 
wave  reflection  coefficients  over  a  range  of  about  15-65  kHz.  Such  data  provide  a  means  for  developing 
phenomenological  models  of  the  lower  ionosphere,  useful  for  VLF/LF  propagation  prediction 
purposes. 

3.  DATA  DISPLAY  AND  REDUCTION  FORMATS 

Three  formats  have  been  chosen  for  representing  the  data  obtained  during  the  experiments,  and 
for  describing  the  VLF/LF  reflection  properties  of  the  lower  ionosphere  that  were  derived  from  it.  A 
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Figure  3.  Basic  VLF/LF  Ionosoundlng  Technique. 
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Figure  4.  Orthogonal  Loop  Antennas  Used  During  Experiment  (Small  Loop 
Used  to  Induce  Calibration  Signals).  Site  Pictured  was  at  Blumenau.  Brazil. 
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Figure  5.  Example  of  Observed  Waveforms. 
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Figure  6.  Fourier  Spectrum  for  Groundwave  Pulse  Shown  in  Figure  5. 
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brief  description  and  a  few  samples  of  the  data  formats  are  given  below;  following  that,  data  obtained 
at  each  of  the  nine  receiving  locations  are  presented. 


3.1.  Three-Dimensional  Data  Displays 

Figure  7  shows  a  three-dimensional  presentation  of  a  sample  of  data  to  illustrate  the  manner  in 
which  data  from  each  receiving  location  were  acquired  and  initially  displayed.  The  presentation 
consists  of  waveforms  stacked  one-behind-the-other,  with  time  progressing  linearly  from  bottom-to- 
top,  for  (in  this  example)  a  three-day  period.  Although  the  raw  data  were  recorded  every  2.8  minutes, 
for  display  purposes  each  waveform  in  the  figure  represents  a  15  minute-average  of  approximately 
325.000  pulses.  The  horizontal  scale  is  linear  in  time,  measured  from  the  start  of  the  groundwaves. 
The  groundwaves  remain  fixed  In  time,  while  the  skywaves  move  back  and  forth  in  accordance  with 
the  solar  zenith  angle  (and  hence,  In  accordance  with  the  varying  effective  height  of  the  lower 
ionosphere).  The  waveforms  represent  the  averaged  wave  amplitudes  as  a  function  of  time,  as  sensed 
by  the  receiving-loop  antennas.  The  data  in  Figure  7  correspond  to  the  normal  component  of  the 
ionospherically  reflected  waves.  Similar  displays  are  given  later  showing  converted  components  of 
the  reflected  waves  received  at  a  number  of  locations  during  the  experimental  campaign. 


3.2  Reflection  Heights  And  Coefficients  As  A  Function  Of  Time 

Data,  such  as  that  illustrated  In  Figure  7.  provide  the  bases  for  deriving  effective  heights  of 
reflection  and  plane-wave  reflection  coefficients  of  the  Ionosphere,  as  outlined  in  Section  2  of  this 
report.  Figure  8  gives  examples  of  reflection  coefficients,  plotted  as  a  function  of  time  for  15,  30.  45, 
and  60  kHz,  as  derived  from  the  sample  of  data  given  in  Figure  7.  In  doing  this,  the  received  skywave 
data  were  averaged  over  2-hr  intervals.  Such  data  presentations  are  given  to  show  longer-  term  trends 
In  the  data  (days  or  weeks),  rather  than  short-term  characteristics  (minutes  or  hours).  Also  shown  In 
Figure  8  Is  a  plot  of  the  effective  height  of  the  lower  ionosphere  at  30.5  kHz.  as  derived  from  the  same 
data. 


3.3  Reflection  Heights  And  Coefficients  As  A  Function  Of  Frequency 

By  fixing  a  time  (for  example,  local  noon) ,  the  pulse  reflection  data  can  be  used  to  derive  effective 
heights  of  reflection  and  plane-wave  reflection  coefficients  as  a  function  of  frequency.  For  example. 
Figure  9  Illustrates  the  results  of  such  an  analysis,  derived  from  the  first  day's  data  at  local  midnight 
and  local  noon,  as  shown  In  Figure  7.  The  data  provide  reflection  properties  of  the  lower  ionosphere 
over  a  wide  frequency  range,  from  about  15-65  kHz.  The  resolution  of  the  system  Is  such  that  a  change 
in  the  effective  height  of  the  Ionosphere  by  as  little  as  0.5  km  can  be  discerned.  In  addition,  reflection 
coefficients  as  small  as  about  0.02  In  amplitude  can  be  determined.  Data  such  as  that  shown  in  Figure 
9  were  derived  from  15-min  averages  of  the  raw  data. 
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Figure  7.  Example  of  3-Dimensional  Display  of  Groundwave 
and  Normal  Skywave  Pulses  Received  at  Canasvieiras  Stacked 
Front  to  Back  with  Advancing  Time. 
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Figure  8.  Example  of  Group  Heights  of  Reflection  and  Normal  Reflection  Coefficients 
Derived  from  Data  Presented  in  Figure  7.  Group  Heights  Calculated  from  Data  at  30.5  kHz 
and  Normal  Reflection  Coefficients  Calculated  at  15,  30,  45  and  60  kHz. 
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b.  LOCAL  NOON 


Figure  9.  Example  of  Group  Heights  of  Reflection  and  Normal  Reflection  Coefficients 
vs  Frequency  at  Local  Midnight  (a)  and  Local  Noon  (b)  Derived  from  Data  of  16  May 
(Day  of  Year  137)  1980. 


11 


4.  DATA  FROM  EACH  OF  THE  RECEIVING  LOCATIONS 


In  this  section  the  data  obtained  during  the  experiment  in  Brazil  are  displayed.  The  data  were 
taken  during  the  late  fall  and  early  winter,  when  the  day-to-day  variations  in  the  solar  zenith  angles, 
as  a  function  of  time  of  day,  were  not  large  and  the  periods  of  sunlight  and  darkness  were  about  equal, 
as  shown  in  Figure  10.  In  the  presentations  that  follow,  the  displays  parallel  those  described  in 
Section  3;  specifically:  three-dimensional  displays  of  the  data  are  first  shown;  these  are  followed  by 
plots  of  effective  heights  of  reflection  (at  30.5  kHz)  and  plane-wave  reflection  coefficients  (at  15,  30.  45, 
and  60  kHz),  as  a  function  of  time;  finally,  reflection  heights  and  reflection  coefficients  are  given,  as  a 
function  of  frequency,  for  local  midnight  and  local  noon  for  the  days  data  were  obtained  at  each 
receiving  site.  Table  1  gives  the  geomagnetic  coordinates,  the  great-circle  propagation  distances,  and 
the  magnetic  azimuths  from  the  transmitter  to  each  receiving  location.  The  specific  periods  for  which 
data  were  acquired  at  each  of  the  receiving  locations  are  shown  on  charts  immediately  preceding  the 
data  displays  for  that  location.  Those  charts  also  indicate  if  only  the  normal  components  of  the 
reflected  skywaves  were  received,  or  whether  the  converted  components  were  received  as  well. 


Table  1.  Locations  and  Propogation  Path  Parameters 
of  the  Transmitting  and  Receiving  Locations  for  the  Experiment. 


SITE 

COORDINATES 

DISTANCE 

KM. 

MAGNETIC 

AZIMUTH 
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Figure  10.  Solar  Zenith  Angle  Variations  for  the  Transmitter  at  Paula 
Freitas  During  the  Measurements  Between  8  May  and  27  June,  1980. 
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PROPAGATION  PATH: 

DISTANCE; 

AZIMUTH: 

DATA: 

MEASUREMENT  PERIOD: 


DATA 

PAULA  FREITAS  TO  DR.  PEDRINHO  (PF-PED) 

156  KM 

123° 

NORMAL  AND  CONVERTED  POLARIZATIONS 
17  MAY  -  21  MAY  1980 


Figure  11.  Path  Parameters  for  Paula  Freitas  to  Dr.  Pedrinho. 
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(b)  CONVERTED  SKYWAVE 


Figure  1 1.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection 
Coefficients  at  15,  30,  45  and  60  kHz  from  Data  Received  at  Dr.  Pedrinho. 

a.  Normal  Skywave. 

b.  Converted  Skywave. 
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Figure  1 1.3.  Group  Heights  of  Reflection  and  Reflection  Coefficients 
vs  Frequency  Derived  from  Data  Received  at  Dr.  Pedrinho. 

b.  Local  Noon  (1515  UT) 
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Figure  1 1 .4.  Averaged  Values  of  Group  Heights  of  Reflection  and  Reflection 
Coefficients  vs  Frequency  Derived  from  Data  Received  at  Dr.  Pedrinho. 

a.  Group  Heights  of  Reflection. 

b.  Reflection  Coefficients. 
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DATA 


PROPAGATION  PATH:  PAULA  FREITAS  TO  TIMBO  (PF-TIM) 
DISTANCE;  180  KM 
AZIMUTH:  124° 

DATA:  NORMAL  AND  CONVERTED  POLARIZATIONS 

MEASUREMENT  PERIOD:  13  MAY  -  17  MAY  1980 


Figure  12.  Path  Parameters  for  Paula  Freitas  to  Timbo. 
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Figure  12.1.  Three-Dimensional  Display  of  Data  Received  at  Timbo 


a.  Normal  Skywave. 
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Figure  12.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection 
Coefficients  at  15,  30,  45  and  60  kHz  from  Data  Received  at  Timbo. 

a.  Normal  Skywave. 

b.  Converted  Skywave. 
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Figure  12.3.  Group  Heights  of  Reflection  and  Reflection  Coefficients 
vs  Frequency  Derived  from  Data  Received  at  Timbo. 

b.  Local  Noon  (1515  UT). 
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Figure  12.4.  Averaged  Values  of  Group  Heights  of  Reflection  ami  Reflection 
Coefficients  vs  Frequency  Derived  from  Data  Received  at  Timbo. 
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DATA 


PROPAGATION  PATH:  PAULA  FREITAS  TO  BLUMENAU  (PF-BLU) 
DISTANCE;  197  KM 
AZIMUTH:  122° 

DATA:  NORMAL  AND  CONVERTED  POLARIZATIONS 

MEASUREMENT  PERIOD:  21  MAY  -  4  JUNE  1980 


Figure  13.  Path  Parameters  for  Paula  Freitas  to  Blumenau. 
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a.  Normal  Skywavo. 
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Figure  13.1.  Three-Dimensional  Display  of  Data  Received  at  Blumenau 


b.  Converted  Skywave 
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Figure  13.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection 
Coefficients  at  15,  30,  45  and  60  kHz  from  Data  Received  at  Blumenau. 

a.  Normal  Skywave. 

b.  Converted  Skywave. 
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Figure  13.3.  Group  Heights  of  Reflection  and  Reflection  Coefficients 
vs  Frequency  Derived  from  Data  Received  at  Blumenau. 

a.  Local  Midnight  (0315  UT). 
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DATA 


DISTANCE;  215  KM 
AZIMUTH:  123° 

DATA:  NORMAL  POLARIZATION  ONLY 

MEASUREMENT  PERIOD:  8  MAY  -  15  MAY  1980 

1  JUN  -  27  JUN  1980 


Figure  14.  Path  Parameters  for  Paula  Freitas  to  Gaspar. 


39 


Figure  14.1.  Three-Dimensional  Display  of  Normal  Skywave  Data  Received  at  Gaspar.  (Coni.) 
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Figure  14.1.  Three-Dimensional  Display  of  Normal  Skywave  Data  Received  at  Gaspar.  (Cont.) 
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Figure  14.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection  Coefficients 
at  15,  30,  45  and  60  kHz  Derived  from  Normal  Skywave  Data  Received  at  Gaspar. 
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Figure  14.4.  Averaged  Values  of  Group  Heights  of  Reflection  and  Reflection 
Coefficients  vs  Frequency  Derived  from  Data  Received  at  Gaspar. 

a.  Group  Heights  of  Reflection. 

b.  Reflection  Coefficients. 
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Figure  15.  Path  Parameters  for  Paula  Freitas  to  Camboriu. 
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Figure  15.1.  Three-Dimensional  Display  of  Normal  Skywave  Data  Received  at  Camboriu.  (Cont.) 
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Figure  15.1.  Three-Dimensional  Display  of  Normal  Skywave  Data  Received  at  Camboriu.  (Cont 
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Figure  15.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection  Coefficients 
at  15.  30,  45  and  60  kHz  Derived  from  Normal  Skywave  Data  Received  at  Camboriu. 
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Figure  15.3.  Group  Heights  of  Reflection  and  Reflection  Coeflicients 
vs  Frequency  Derived  from  Data  Received  at  Camboriu. 


a.  Local  Midnight  (0315  UT). 

b.  Local  Noon  (1515  UT). 


Figure  15.4.  Averaged  Values  of  Group  Heights  of  Reflection  and  Reflection 
Coefficients  vs  Frequency  Derived  from  Data  Received  at  Camboriu. 

a.  Group  Heights  of  Reflection. 

b.  Reflection  Coefficients. 
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DATA 


PROPAGATION  PATH:  PAULA  FREITAS  TO  CANASVIEIRAS  (PF-CAN) 
DISTANCE;  280  KM 
AZIMUTH:  132° 

DATA:  NORMAL  POLARIZATION  ONLY 

MEASUREMENT  PERIOD:  15  MAY  -  19  MAY  1980 


Figure  16.  Path  Parameters  for  Paula  Freitas  to  Canasvieiras. 
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Figure  16.1.  Three-Dimensional  Display  of  Normal  Skywave  Data  Received  at  Canasvieiras. 
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Figure  16.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection  Coefficients  at  15,  30.  45 
and  60  kHz  Derived  from  Normal  Skywave  Data  Received  at  Canasvieiras. 
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Figure  16.3.  Group  Heights  of  Reflection  and  Reflection  Coefficients 
vs  Frequency  Derived  from  Data  Received  at  Canasvieiras. 


a.  Local  Midnight  (0315  UT). 

b.  Local  Noon  (1515  UT). 
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Figure  16.4.  Averaged  Values  of  Group  Heights  of  Reflection  and  Reflection 
Coefficients  vs  Frequency  Derived  from  Data  Received  at  Canasvieiras. 

a.  Group  Heights  of  Reflection. 

b.  Reflection  Coefficients. 
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Figure  17.  Path  Parameters  for  Paula  Freitas  to  Figueira. 
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Figure  17.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection  Coefficients  at 
15,  30,  45  and  60  kHz  Derived  from  Normal  Skywave  Data  Received  at  Figucira. 
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Figure  17.3.  Group  Heights  of  Reflection  and  Reflection  Coefficients 
vs  Frequency  Derived  from  Data  Received  at  Figueira. 

a.  Local  Midnight  (0315  UT). 

b.  Local  Noon  (1515  UT). 
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Figure  17.4.  Averaged  Values  of  Group  Heights  of  Reflection  and  Reflection 
Coefficients  vs  Frequency  Derived  from  Data  Received  at  Figueira. 

a.  Group  Heights  of  Reflection. 

b.  Reflection  Coefficients. 
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DATA 

PROPAGATION  PATH:  PAULA  FREITAS  TO  PINHALZINHO  (PF-PIN) 


DISTANCE:  213  KM 

AZIMUTH:  262° 

DATA:  NORMAL  POLARIZATION  ONLY 

MEASUREMENT  PERIOD:  29  MAY  -  1  JUN  1980 


Figure  18.  Path  Parameters  for  Paula  Freitas  to  Pinhalzinho. 
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Figure  18.1.  Three-Dimensional  Display  of  Normal  Skywave  Data  Received  at  Pinhalzinho. 
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Figure  18.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection  Coefficients  at 
15,  30,  45  and  60  kHz  Derived  from  Normal  Skywave  Data  Received  at  Pinhalzinho. 
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Figure  18.4.  Averaged  Values  of  Group  Heights  of  Reflection  and  Reflection 
Coefficients  vs  Frequency  Derived  from  Data  Received  at  Pinhalzinho. 


a.  Group  Heights  of  Reflection. 

b.  Reflection  Coefficients. 
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PROPAGATION  PATH:  PAULA  FREITAS  TO  IRAI  (PF-IRA) 
DISTANCE;  250  KM 
AZIMUTH:  256° 

DATA:  NORMAL  POLARIZATION  ONLY 

MEASUREMENT  PERIOD:  26  MAY  -  29  MAY  1980 


Figure  19.  Path  Parameters  for  Paula  Freitas  to  Irat. 
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Figure  19.1.  Three-Dimensional  Display  of  Normal  Skywave  Data  Received  at  Irai. 
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Figure  19.2.  Group  Heights  of  Reflection  for  30.5  kHz  and  Reflection  Coefficients 
at  1 5.  30,  45  and  60  kHz  Derived  from  Normal  Skywave  Data  Received  at  Irai. 
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Figure  19.4.  Averaged  Values  of  Group  Heights  of  Reflection  and  Reflection 
Coefficients  vs  Frequency  Derived  from  Data  Recieved  at  Irai. 


a.  Group  Heights  of  Rpflection. 

b.  Reflection  Coefficients. 


5.  SELECTED  FEATURES  OF  THE  REFLECTION  DATA 


In  this  section  data  are  shown  to  illustrate  variations  in  the  VLF/LF  reflection  properties  of  the 
ionosphere  with  such  factors  as  time-of-day.  direction  of  propagation,  range  and  skywave 
polarization.  In  addition,  a  brief  description  is  given  of  weak  daytime  pulses,  seen  from  time-to-time, 
which  were  reflected  from  altitudes  somewhat  lower  than  those  normally  observed. 

Due  to  the  limited  amount  of  data  obtained  at  many  of  the  receiving  locations,  the  variations 
described  below  should  be  considered  to  be  indicative  of  trends,  rather  than  more-fully  established 
characterizations.  These  trends,  however,  are  consistent  generally  with  those  observed  in  earlier 
experimental  programs,  conducted  at  both  mid-  and  low-geomagnetic  latitudes  2' 1 . 

5.1  Day/Night  Reflection  Heights 

A  summary  of  noon  and  midnight  reflection  heights,  derived  from  the  data  obtained  at  all  the 
receiving  locations,  is  given  in  Figure  20.  As  shown,  the  effective  heights  of  reflection  varied  little 
over  the  frequency  range  from  15  kHz  to  60  kHz.  At  noon  the  average  height  of  reflection  was  69  km;  at 
midnight,  it  was  86  km.  Throughout  the  15-60  kHz  range,  the  spreads  in  the  heights  were  relatively 
small,  usually  well  less  than  2.5  km  and  1.5  km  for  the  noon  and  midnight  data,  respectively. 


5.2  Day/Night  Variations  In  The  Reflection  Coefficients 

Normal  (||R||)  reflection  coefficients,  derived  from  noon  and  midnight  data  obtained  at  a  number 
of  receiving  locations  to  the  east  and  west  of  the  transmitter,  are  given  in  Figure  21.  In  all  cases  the 
reflection  coefficients  at  midnight  were  substantially  larger  than  those  at  noon,  ranging  from  a  factor 
of  about  2  for  the  lowest  frequencies  observed  (15  kHz)  to  an  order  of  magnitude  or  more  at  much 
higher  frequencies  (40-60  kHz).  At  midnight  and  at  noon  the  magnitudes  of  the  reflection  coefficients 
generally  decreased  as  the  frequency  increased.  At  midnight,  however,  the  variations  were  relatively 
small,  usually  being  only  about  a  factor  or  two  over  the  15-60  kHz  band.  In  contrast,  at  noon  the 
magnitudes  of  the  reflection  coefficients  varied  by  an  order  of  magnitude  or  more  over  the  same 
frequency  range.  These  results  suggest  that  at  midnight  the  reflections  effectively  came  from  a 
sharply  bounded  layer,  while  at  noon  there  was  weak  ionization,  below  the  altitudes  at  which  the 
waves  were  reflected,  that  absorbed  appreciable  energy  from  the  up-  and  down-going  waves.  The  fact 
that  the  effective  heights  of  reflection  varied  so  little  over  the  15-60  kHz  band  (see  Section  5.1)  also 
supports  such  a  characterization  of  the  difference  between  the  midnight  and  noon  reflection 
processes. 

5.3  Azimuthal  Variations  In  The  Reflection  Coefficients 

Comparison  of  the  east-west  (magnetic)  data  in  Figure  21  reveals  that  the  reflection  coefficients 
at  noon  were  appreciably  larger  for  propagation  to  the  east.  This  directional  effect,  which  is  due  to  the 
influence  of  the  geomagnetic  field,  is  in  accordance  with  earlier  experimental  observations2  and 
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GROUP  MIRROR  HEIGHT  -  km 


Figure  20.  Summary  of  Noon  and  Midnight  Group  Heights  of  Reflection 
Derived  from  Data  Obtained  at  all  the  Receiving  Locations. 
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Figure  21.  Normal  Reflection  Coefficients  for  Midnight  and  Noon  Derived  from 
Data  at  Selected  Locations  to  the  East  and  to  the  West  of  the  Transmitter. 

a.  Approximately  200  km  Propagation  Paths  (Paula  Freitas  to  Blumenau/Figueira). 

b.  Approximately  215  km  Propagation  Paths  (Paula  Freitas  to  Gaspar/Pinhalzinho) 

c.  Approximately  250  km  Propagation  Paths  (Paula  Freitas  to  Camboriu/lrai). 
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theoretical  expectations5.  The  midnight  data  exhibit  little  east-west  difference,  even  though  the 
reflections  came  from  much  higher  altitudes  where  such  geomagnetic  effects  would  normally  be 
expected  to  be  stronger,  because  of  the  reduced  collisional  forces  associated  with  the  lower  neutral 
densities  at  those  heights.  The  absence  of  appreciable  geomagnetic  effects  in  the  midnight  data 
suggests  again,  as  described  in  Sections  5. 1  and  5.2.  that  the  nighttime  reflection  process  can  be 
characterized  as  having  been  more  sharply  bounded  than  the  daytime  process. 


5.4  Variations  In  The  Reflection  Coefficients  With  Range 

Figures  22  and  23  summarize  the  variations  in  the  magnitudes  of  the  normal  reflection 
coefficients  with  range  {or,  in  effect,  incidence  angle)  that  were  derived  from  the  pulse  reflection  data. 
In  Figure  22  data  are  given  at  four  frequencies  over  (he  15-60  kHz  band  for  six  receiving  locations, 
ranging  from  about  156  to  about  280  km  to  the  east  of  the  transmitter.  This  corresponds,  effectively, 
to  a  measure  of  the  reflectivity  of  the  ionosphere  over  a  range  of  Incidence  angles  from  about  48-64 
degrees  in  the  day.  or,  about  42  48  degrees  at  night.  The  data  show  that  the  magnitudes  of  the 
reflection  coefficients  generally  increase  with  increasing  range,  or  more  grazing  incidence  angle,  with 
the  trend  being  especially  pronounced  in  the  daytime. 

Figure  23  shows  20  kHz  reflection  coefficients  derived  from  data  obtained  at  locations  to  the  east 
and  west  of  the  transmitter.  To  the  east,  the  reflection  coefficient  Increases  monotonically  as  the 
length  of  the  propagation  path  increases.  To  the  west,  however,  there  is  an  apparent  dip  in  the 
coefficient  at  a  range  of  about  215  km  (corresponding  to  an  incidence  angle  of  57  degrees),  which  may 
be  indicative  ol  a  quasi-Brewster’s  angle  effect  in  the  ionosphere.  These  data  agree  almost  exactly  with 
similar  data  obtained  near  the  geomagnetic  equator  in  Brazil  in  an  earlier  experimental  program1 
(Rasmussen  et  al). 


5.5  Comparison  Of  Normal  And  Converted  Reflection  Coefficients 

Only  a  small  amount  of  data  pertaining  to  the  polarization  rotation  effects  produced  in  the 
ionosphere  by  the  geomagnetic  field  were  acquired  during  the  experimental  campaign.  The  results  of 
those  observations  arc  summarized  in  Figure  24,  which  gives  normal  (i|Rn)  and  converted  (iiRi) 
reflection  coefficients  derived  from  noon  and  midnight  data  obtained  at  three  locations  to  the  east  of 
the  transmitter.  As  shown,  the  normal  reflection  coefficients  were  larger  than  the  converted 
coefficients,  with  the  difference  being  especially  great  during  the  daytime  for  frequencies  below  about 
30  kHz  In  general,  the  polarization  conversion  at  night  was  substantially  greater  than  in  the 
daytime,  as  expected  from  considerat  of  the  relative  effects  of  electron-neutral  collisions  in  the 
reflection  pro<  ess.  Theoretical  considerations  also  indicate,  generally,  that  while  the  normal 
reflection  coefficients  should  inr lease  as  the  dist-anc'>s  from  the  transmitter  become  greater  and 
greater  (as  drsi  (Ted  m  Section  f>.  11.  the  . inversion  cneflkier's  should  decrease.  Unfortunately,  the 


5.  Crombie,  D  U  (1358)  Differences  between  V.  \V  and  WE  propagation  of  VLF  signals  over  long 
distances.  •  /.  At mos.  Ten.  Phns  12-  1 10 
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Figure  22.  Variation  in  Normal  Reflection  Coefficients  with  Range 
and  Angle  of  Incidence  for  Locations  East  of  the  Transmitter. 
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Figure  23.  Normal  Reflection  Coefficients  for  20  kHz  Derived  from  Midnight  and  Noon 
Data  Obtained  at  Receiving  Locations  to  the  East  and  West  of  the  Transmitter. 
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Figure  24.  Normal  and  Converted  Reflection  Coefficients  Derived  from 
Midnight  and  Noon  Data  from  Three  Locations  East  of  the  Transmitter. 
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polarization  rotation  data  that  was  acquired  during  the  experimental  campaign  was  too  sparse  to  test 
this  theoretical  expectation  adequately. 

5.6  Evidence  Of  Weak  Reflections  From  Below  The  Normal  D-Region 

The  data  acquired  during  sunrise  periods  were  inspected  for  evidence  of  reflections  coming  from 
altitudes  below  the  normal  D-region,  such  as  those  observed  earlier  in  similar  mid-latitude 
experiments  6,  7  which  were  attributed  to  a  'C-layer'  of  ionization  produced  by  the  combined  action  of 
cosmic  rays  and  solar  illumination.  A  total  of  57  sunrise  recordings  were  made  during  the  22  days 
comprising  the  experimental  campaign.  Only  rarely  were  recordings  made  during  any  given  sunrise  at 
more  than  two  locations.  Of  the  57  sunrise  recordings,  approximately  40  percent  showed  evidence  of 
reflections  below  the  normal  D-region.  In  most  cases,  however,  the  reflections  were  very  weak  and 
partially  mixed  with  the  D-region  reflections;  that  is,  only  the  leading  edge  of  the  received  pulses  could 
be  resolved  from  the  later  arriving,  D-region.  reflections.  In  addition,  there  were  sunrises  for  which 
the  low  altitude  reflections  were  observed  at  one  receiving  location,  but  not  another. 

Figure  25  shows  data  acquired  over  a  four-day  period  at  Camboriu,  246  km  to  the  east  of  the 
transmitter.  Examples  of  weak  reflections,  believed  to  have  come  from  a  C-layer,  are  outlined  by 
circles  in  the  Figure.  Such  reflections  are  most  easily  seen  in  the  data  of  May  17,  where  the  leading 
edges  of  the  daytime  skywaves  were  nearly  stationary,  and  distinct  from  later  arriving  (D-region) 
reflections  that  moved  in  accordance  with  the  solar  zenith  angle.  However,  on  other  days  shown  in 
the  Figure,  such  distinct  reflections  could  only  be  seen  for  very  short  periods  around  sunrise  and 
sunset. 

Very  clear  evidence  of  the  reflections  occurred  during  sunrise  on  May  17,  in  observations  made 
simultaneously  at  Timbo,  Camboriu  and  Canasvieiras,  at  distances  of  180,  246,  and  280  km, 
respectively,  to  the  east  of  the  transmitter.  The  pulse  reflection  data  recorded  at  those  sites  at  1015  UT, 
when  the  solar  zenith  angle  was  about  86  degrees,  is  given  in  Figure  26.  In  the  three  traces  of  Figure 
26a,  two  skywave  pulses  are  seen,  which  are  clearly  separated  from  an  earlier  arriving  groundwave 
pulse.  The  first  skywave  (labeled  C  in  each  tract)  is  presumed  to  have  been  reflected  from  a  C-layer  in 
the  lower  ionosphere,  whereas  the  second,  much  stronger,  skywave  is  presumed  to  have  been  reflected 
from  D-region  ionization  produced  by  Lyman-a  radiation  from  the  sun.  Based  on  the  arrival  times  of 
the  early  skywave  reflections  and  the  path  lengths  over  which  the  observations  were  made  it  Is 
estimated  that  the  effective  height  of  the  C-layer  that  produced  the  reflections  was  67  km. 

Figure  26b  shows  effective  plane-wave  reflection  coefficients  derived  from  the  C-layer 
reflections.  The  coefficients  provide  estimates  of  the  reflection  properties  of  the  C-layer  as  a  function 
of  frequency  and  incidence  angle,  which  can  then  be  used  to  infer  a  conducting-slab  model  of  the  layer. 


6.  Rasmussen,  J.E.,  Kossey,  P.A.,  and  Lewis,  EA  (1980)  Evidence  of  an  ionospheric  reflecting  layer 

below  the  classical  D  region,  J.  Geophys.  Res.  85:  3037-3044. 

7.  Bain,  W.C.  and  Kossey,  P.A  (1987)  Characteristics  of  a  reflecting  layer  below  the  classical  D 

region"  J.  Geophys.  Res.  92:  12,443-12,444. 
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Figure  26.  Data  Recorded  on  17  May  (Day  138)  1980  at  Timbo.  Camboriu  and  Canasvieiras. 

a.  Normal  Skywave  Data  Showing  Weak  Reflections  below  the  Normal  D-Region  (C). 

b.  Normal  Reflection  Coefficients  Derived  from  the  "C”  Protions  of  the  Received  Skywaves. 
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As  described  by  Kossey  and  Lewis  8,  an  estimate  of  the  thickness  of  the  slab  can  be  obtained  from 
the  frequencies  at  which  the  relative  nulls  occur  in  the  reflection  coefficients.  The  magnitudes  of  the 
reflection  coefficients  are  then  used  in  determining  an  effective  conductivity  of  the  slab.  From  the 
data  of  Figure  26b,  it  was  estimated  that  a  10-km  thick  slab,  having  a  uniform  conductivity  of  10'7 
Siemens/m,  could  describe  the  ionization  (phenomenologically)  that  produced  the  C-layer  reflections. 
This  model  contrasts  with  a  more  strongly  ionized.  6-km  thick  slab,  located  at  60  km  altitude  that  was 
inferred  from  earlier  pulse  reflection  measurements  at  a  mid-latitude  location  in  the  United  States  7 . 

If  the  low  altitude  reflections  that  were  observed  in  Brazil  did  indeed  come  from  a  layer 
produce/ d  by  cosmic  rays,  it  is  not  too  surprising  that  its  features  might  differ  appreciably  from  those 
observed  at  higher  latitudes.  This  follows,  qualitatively  at  least,  from  consideration  of  the 
geomagnetic  field  conditions.  In  Brazil,  the  measurements  were  made  where  the  geomagnetic  field  had 
a  dip-angle  of  only  about  15  degrees;  that  is,  the  field  was  almost  horizontal.  It  is  surmised  that  'he 
nearly  horizontal  geomagnetic  field  produces  a  shielding  effect  to  charged  particles,  thus  reducing  the 
ability  of  the  cosmic  rays  to  penetrate  deeply  into  the  lower  atmosphere.  At  higher  latitudes,  where 
the  geomagnetic  field  has  a  much  larger  dip-angle,  such  shielding  effects  are  expected  to  be  veiy  much 
weaker.  For  example,  at  the  mid-latitude  location  where  the  pulse  reflection  data  were  used  to  infer  a 
relatively  strongly  ionized,  6-km  thick  layer,  at  60  km  altitude,  the  dip-angle  of  the  geomagnetic  field 
was  almost  70  degrees. 

Too  little  appropriate  data  were  acquired  during  the  experiments  in  Brazil  to  attempt  to 
characterize  the  nature  of  the  C-layer  at  low  latitudes.  The  limited  data  that  were  obtained,  however, 
provide  useful  information  on  how  future  experiments  might  be  designed  to  better  exploit  the  pulse 
reflection  technique  for  investigating  the  C-layer.  For  example,  the  data  of  Figure  26  illustrate  that 
the  lengths  of  the  propagation  paths  in  such  experiments  should  be  carefully  chosen.  If  the  path  is  too 
short,  (that  is,  the  incidence  angle  of  the  pulses  on  the  layer  is  too  steep)  the  reflections  may  be  too 
weak  to  be  clearly  discerned.  For  example,  the  C-layer  reflection  observed  in  Timbo,  180  km  from  the 
transmitter,  was  almost  too  weak  to  be  measured;  while  in  Canasvieiras,  at  a  distance  of  280  km,  it 
was  relatively  strong,  and  easily  measured.  On  the  other  hand,  at  Canasvieiras,  there  was  very  little 
separation  between  the  groundwave,  C-layer  and  D-region  pulse  reflections,  owing  to  the  relatively 
long  propagation  path.  At  the  intermediate  distance  of  246  km,  at  Camboriu,  however,  there  was 
ample  separation  between  the  observed  pulses,  and  the  amplitude  of  the  C-layer  reflection  was 
relatively  strong. 


5.7  Unusual  Moving  High  Altitude  Nighttime  Reflections 

Often  at  night  unusual,  moving,  reflections  were  observed  coming  from  regions  well  above  the 
altitudes  normally  associated  with  the  VLF/LF  pulse  reflections.  The  reflections  tended  to  begin  at  a 
high  altitude  around  0200  UT  and,  in  most  cases,  drifted  downward  over  a  period  of  6-8  hr  before 
disappearing  during  the  sunrise  transition  period.  Sometimes,  however,  the  reflections  began  to  drift 
upward  again,  about  2  hr  before  sunrise.  These  moving  reflections  can  be  easily  seen  in  many  of  the 


8.  Kossey,  P.A.  and  Lewis,  E.A.  (1981)  The  Reflection  Properties  of  Conducting  Slabs,  RADC-TR-80- 
371,  ADA098940. 
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quasi  3-dimensional  displays  of  reflection  data  given  in  this  report.  A  five-day  sample  of  such  data  is 
given  in  Figure  27.  It  is  hoped  that  future  research  will  provide  the  necessary  analyses  to  characterize 
these  unusual  nighttime  reflections  and  the  ionization  processes  that  produce  them. 


6.  ELECTRON  DENSITY  MODELS  OF  THE  LOWER  IONOSPHERE 

Phenomenological  models  of  the  lower  ionosphere  can  be  derived  from  steep-incidence  VLF/LF 
reflection  data  using  mathematical  inversion  techniques  such  as  those  summarized  by  Sechrist9. 
Under  normal  ionospheric  conditions  this  is  not  an  easy  task,  since  the  polarization  rotation  effects 
of  the  geomagnetic  field  cannot  usually  be  ignored.  In  addition,  there  is  a  question  as  to  the 
uniqueness  of  models  derived  from  reflection  data.  Nevertheless,  since  the  VLF/LF  reflection 
properties  are  very  sensitive  to  electron  densities  below  about  90  km  at  night,  and  below  75  km  during 
the  day,  such  data  provide  an  important  means  for  developing  and  validating  models  of  the  lower 
ionosphere. 

Figure  28  shows  electron-density  models  of  the  low  latitude  ionosphere  in  Brazil,  derived  from 
the  pulse  reflection  data.  The  models  were  obtained  using  a  trial-and-error  approach,  which  employed 
full-wave  and  iterative  computational  techniques.  In  doing  this,  data  from  all  nine  of  the  receiving 
sites,  were  included  so  that  the  derived  models  were  made  to  produce  computational  results  consistent 
with  a  large  number  and  variety  of  experimental  data. 

The  technique  involved  the  choice  and  refinement  of  trial  electron  density  models  of  the  lower 
ionosphere  that  served  as  inputs  to  a  full-wave  computational  program,  used  to  calculate  theoretical 
values  of  VLF/LF  reflection  coefficients  and  effective  heights  of  reflection,  across  a  wide  frequency 
band.  The  calculated  values  were  then  compared  to  those  derived  from  the  experimental  data,  and  the 
differences  between  them  were  used  as  the  basis  for  adjusting  the  trial  electron  density  profiles. 

Finally,  for  each  electron  density  model  a  root-mean-square  (RMS)  error  was  determined  from  the 
theoretical  (computer)  and  experimental  values.  As  such,  the  RMS  errors  provided  a  means  for 
comparing  the  various  electron  density  models  under  consideration. 

Based  on  the  strong  tendency  for  the  experimental  values  of  the  reflection  coefficients  to 
decrease  monotonically  with  increasing  frequency,  trial  electron-density  profiles  were  chosen  that 
varied  exponentially  with  altitude.  Figure  28  shows  the  noon  and  midnight  electron-density  models, 
having  the  smallest  RMS  errors,  that  were  obtained  from  the  numerous  exponential  models  that  were 
tested.  In  the  notation  of  Wait  and  Spies  10,  which  is  often  used  to  characterize  such  models,  the  noon 
model  is  described  by  a  (3  =  0.43  km-1  and  a  h  =  70  km,  and  the  midnight  model  by  (3  =  0.8  km-1  and  h  = 
83  km.  For  the  noon  model,  the  average  RMS  error  between  the  calculated  and  experimental  values  for 
the  magnitudes  of  the  normal  reflection  coefficients  was  1.4  dB  over  the  15-50  kHz  band  and  averaged 
over  the  nine  receiving  locations  in  Brazil.  The  corresponding  RMS  error  in  the  eilective  heights  of 
reflection  was  1.2  km.  For  the  midnight  model,  the  RMS  errors  were  1.8  dB  and  1.0  km,  respectively. 


9.  Sechrist,  C.F.  (1974)  Comparison  of  techniques  for  measui  ement  of  D  region  electron  densities. 
Radio  Sci  9:  137-149. 

10.  Wait,  J.R.  and  Spies,  K.P.  (1964)  Characterislics  of  the  earth-ionosphere  waveguide  for  VLF  radio 
waves,  Nat.  Bur.  Stds.,  Tech.  Note  300. 
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Figure  27.  Data  Obtained  at  Camboriu  Showing  Moving  Reflections  (circled  areas) 
During  the  Nighttime  coming  from  Relatively  High  Altitudes  above  the  D-Region. 


Figure  28.  Electron  Density  Profiles  of  Low  Latitude  Ionosphere  Derived  from 
the  VLF/LF  Reflection  Data  Gathered  During  this  Series  of  Measurements. 
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determined  from  data  over  an  even  wider  frequency  band  of  15-65  kHz.  The  agreements  between  the 
experimental  and  calculated  values  for  the  two  models  are  very  good,  based  on  earlier  experiences  with 
developing  such  phenomenological  models  of  the  ionosphere  for  example,  Lewis  et  al2  . 

7.  SUMMARY 

Observations  of  VLF/LF  pulse  reflections  from  the  lower  ionosphere,  acquired  at  a  relatively  low 
geomagnetic  latitude  in  Brazil,  have  been  presented.  In  addition,  the  use  of  the  data  to  determine 
effective  heights  of  reflection  and  to  provide  a  measure  of  the  plane-wave  reflection  coefficients  of  the 
ionos’  here,  derived  from  such  data,  were  discussed.  Although  the  models  are  phenomenological,  and 
not  mathematically  unique,  they  should  be  useful  as  inputs  in  longwave  propagation  prediction  codes. 
More  generally,  the  VLF/LF  reflection  data  given  in  this  report  should  be  useful  for  validating  other 
electron  density  models  of  the  lower  ionosphere,  such  as  those  derived  from  long  path  VLF/LF 
radiowave  observations  or  from  theoretical  considerations  of  the  chemistry  of  the  low  latitude 
atmosphere;  that  is,  such  models  should  have  VLF/LF  reflection  properties  that  are  consistent  with 
experimental  data  such  as  the  data  described  in  this  report. 

Other  reflection  data  acquired  during  the  experiments  require  further  study  to  characterize  their 
properties  and  to  identify  the  processes  that  may  be  responsible  for  them.  These  include  weak  daytime 
reflections  from  layers  below  the  classical  D-region,  seen  during  sunrise  periods  and  moving 
reflections  from  altitudes  well  above  those  normally  expected  with  VLF/LF  waves,  observed  from  time 
to  time,  at  night. 
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